Symmetrically coupled higher-order nonlinear 
Schrodinger equations: singularity analysis and 
integrability 

S. Yu. SakovichI and Takayuki Tsuchida| 

f Institute of Physics, National Academy of Sciences, P.O. 72, Minsk, Belarus 

I Department of Physics, Graduate School of Science, University of Tokyo, Hongo 

7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan 

Abstract. The integrability of a system of two symmetrically coupled higher-order 
nonlinear Schrodinger equations with parameter coefficients is tested by means of 
the singularity analysis. It is proven that the system passes the Painleve test for 
integrability only in ten distinct cases, of which two are new. For one of the new 
cases, a Lax pair and a multi-field generalization are obtained; for the other one, the 
equations of the system are uncoupled by a nonlinear transformation. 



f E-mail: sakovichOdragon.bas-net .by 

I E-mail: tsuchidaOmonet . phys . s .u-tokyo . ac. jp 



1. Introduction 



The nonlinear Sclirodinger (NLS) equation, which describes the time evolution of 
slowly varying envelope, is encountered in various branches of physics and known to 
be fundamental. In deriving the NLS equation as the envelope equation, we neglect 
higher-order terms under appropriate physical assumptions. However, due to recent 
developments in optical technology, higher-order corrections to the NLS equation have 
become necessary and important. Kodama and Hasegawa ^ proposed the higher- 
order NLS (HONLS) equation, 

qt = hq^^^ + aqqqx + hq^q^^ + i{sq^x + fq^q) (LI) 

which describes ultra-short pulse propagation in optical fibers, including higher-order 
effects such as higher-order dispersion, self-steeping and delayed Raman response. Here 
the bar denotes the complex conjugation, q = q*. For simplicity, we assume that the 
parameters h, a, b, s, f are real and satisfy the conditions, + b"^ ^ 0, s ^ if h = 0. 
Integrable cases of the HONLS equation ( |1.1|) attract both theoretical and experimental 
interest because they support a variety of exact solutions and the initial-value problem 
is solvable. To date, there are four known integrable cases of the HONLS equation ( |1 . 1| ) 
which are solvable via the inverse scattering method: 

h^O a^O b = f = sa/{3h); (1.2) 

h^O a^O 6 = a/3 f = 2sa/{9h); (1.3) 

h = s^O a^O 6 = 0; (1.4) 

h = s^O a^O 6 = a/2. (1.5) 



The cases ( |1.2|) -(pT5 D of the HONLS equation (|T]T|) are called, respectively, the Hirota 
equation 0], the Sasa-Satsuma equation [§, ||, the Chen-Lee-Liu equation and the 
Kaup-Newell equation 0. The Painleve analysis of the HONLS equation, which has 
been carried out by a number of authors [H, |^, |10|, 0, |12|, strongly indicates that eq. 
( [L.l|) is integrable only in the four cases ( |1.2|) -( pT5D . The study based on prolongation 
structures leads to the same indication [p!3| . 

To describe two pulses copropagating in optical fibers, we need to consider a 



two-component generalization of the single-component propagation equation [T^]. For 
this purpose, we consider the following system of two symmetrically coupled HONLS 
equations: 

qt = hq^^^ + aqqq^ + bq^q^ + crfq^ + dqfr^ + eqrf^ + i(sg^3; + fq^q + gqrf) 

(1.6) 

rt = hr^^^ + arrr^ + br'^r^ + cqqr^ + drqq^, + erqq^ + i(sr^.^. + /r^r + grqq) 

where /i, a, 6, c, e, s, /, g are real parameters and the bar denotes the complex 
conjugation. In fact, we can assume more general conditions q = ±g*, f = ±r* with 
the double signs being arbitrary. This makes no difference in the following analysis. 
We impose the conditions, + b"^ + + d"^ + e"^ ^ 0, + d"^ + + g"^ 0, s if 
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/i = 0, on the parameters so that the system ( |1.6| ) is not uncoupled and includes linear 
dispersion and higher-order nonlinearity. The system ( |1.6| ) is a natural generalization 



of the HONLS equation (LI) to a two- component system, which is invariant under any 



e'^r (/3 : real). 



of the transformations, g r, g — e'"g {a : real), 1 

In this paper, we study the system of coupled HONLS equations ( |L6| ) by means of 
the Painleve analysis. Similar attempts have been performed for coupled NLS equations 
without higher-order terms [|l^, However, on the Painleve analysis of coupled 

HONLS equations, only a few papers of academic significance have appeared Jl^, [T^, 
where the Painleve test has been applied to only integrable cases or a rather restricted 



class of equations, compared with our general form (|L6| ). By using the Painleve test, 
we exhaustively obtain the integrability conditions on the parameters in the coupled 
HONLS equations (|L6| ) for the first time. 

The paper consists of the following. In section H, we perform the singularity analysis 
of (|1.6| ). It is proven that the system ( |L6|) passes the Painleve test for integrability only 
in the following ten distinct cases: 



h^O a^O b=0 c=d 



h^O a^O b=d 
h^O a^O b = e- 



d = g = c 



« f = g = sa/{3h); 
e = a/2 f = sa/{3h); 



h=0 s^O 
h=0 s^O 
h=0 s^O 



h = 
h=0 s^O 



b = d 

a ^ 
a ^ 
s ^ a^O 



e = c = a f = g = sa/ {3h); 

c = d = a/2 f = g = sa/{3h); 

e = a/3 c = 2a/3 f = g = 2sa/{9h); 

b = d = e = c = a g = f] 

b = c = d = e = a/2 g = f] 

b = c = e = d = a g = f] 

b = e = a/2 c = a d = g = f', 

b = d = c = e = a g = 0- 



(L7) 
(L8) 

(1-9) 
1.10) 
1.11) 
1.12) 
1.13) 
1.14) 
1.15) 
1.16) 



In section ^, we show that the system (|1.6| ) is integrable in the Lax sense in all the 
cases (|1.7| )- ([TTT^ ). The integrability of (|1.9|) - (|1.15|) has already been studied in the 
literature. The case ( |1.8|) turns out to be related to ( |1.10|) . For the case (|1.7| ), we 
construct a corresponding 4x4 Lax pair and propose a multi-field generalization. For 
the case (|1.16|) , we obtain a nonlinear transformation, which changes the system ( |1.6| ) 
into two independent Chen-Lee-Liu equations. The last section, section ^, is devoted 
to concluding remarks. 
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2. Singularity analysis 



2.1. Preliminaries 



Let us apply the Painleve test for integrability to the system ( pTq) . It would be better 



to say "the Kovalevskaya test" ||T9|, because we, following the Weiss-Kruskal algorithm 
of the singularity analysis pO|, select only those cases of the tested system, in which 
its general solution admits series expansions of Laurent type. 

With respect to g, g, r, f, which should be considered as mutually independent, the 



r 
f 



(2.1) 



system ( p..6|) is a normal system of four PDEs, of total order twelve if /i 7^ or eight if 
/i = 0, s 7^ 0. A hypersurface t) = is non-characteristic for ( |L6|) if (px 7^ 0, and we 
set 0a; = 1. The substitution of the expansions 

■■■ + g„(t)r+- + --- 

■ ■ ■ + gn(t)0"+^ + • • • 

.■■ + r„(t)0"+^ + --- 

■ ■ + r„(t)0"+^ + ■ ■ ■ 

(the bar does not mean the complex conjugation now) into the system ( |L6|) determines 
the branches, i.e. admissible choices of a, (3, 7, 5 and go? ^o? '"o? "^O; as well as the positions 
n of resonances for those branches. 

We require that the system ( |1.6| ) admits a singular generic branch, where the 
exponents a,/?, 7,5 are integer and at least one of them is negative, the number of 
resonances is equal to the total order of the system, all the resonances but one lie in 
nonnegative integer positions, and the recursion relations for the coefficients of ( |2.1| ) are 
consistent at the resonances. The singular generic branches are studied in sections ^]2| 
and |2]^ for /i 7^ and /i = 0, respectively; other branches are concisely considered in 
section p.4| . Computations are done by means of the Mathematica system ||22|, and we 



omit inessential details for this reason. 



2.2. Singular generic branches: h ^ 

When h 7^ 0, we set /i = 1 without loss of generality. The consideration of 
dominant terms of ( |l.(j| ) leads us to the following two cases to be studied separately: 
a + P = 'J + S = —2 and a + P = —2, 7 + 5 > — 2 (for the reason of symmetry, we omit 
a + /5> -2,7 + 5 = -2). 



2.2.1. Case a + /3 = 7 + 5 = — 2 In this case, we have two possibilities: either two or 
three of the resonances lie in the position n = 0. 

Two resonances at n = 0. If we set a 7^ /5 or 7 7^ S, then the positions of four 
resonances are n = —1,0,0,4, and there are 54 distinct cases of admissible positions 
for other eight resonances. In none of the 54 cases, however, a, (3, 7, 6 are integer 
simultaneously. 
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Therefore we choose a = /3 = 7 = (5 = — 1. Then we find from ( |1.6| ) and ( p.l| ) that 
QoQo = '"o^o = constant 7^ (we set constant = 1 w.l.g.), and that 

a = -6-6-c-rf-e (2.2) 
(n + l)n^{n - 3){n - 4:){n^ -6n-2b-2d + 5){n^ - 6n - 26 - 2e + 5) x 
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6n2 + (5 - 2c/ - 2e)n + 4(c + + e + 3)) = 0. 



:2.3) 



Five resonances he in the positions n = —1,0,0,3,4, and, denoting the positions of 
other seven resonances as ni, n2, . . . , n7, we find from (p.3|) that 



n2 = 6 — rii = 6 — ris rij = Q — — Uq 

d = 1(5-26-6^1 + ^2) e = |(5-26-6n3 + n|) 

h = \{h — Qrii + n\ — 6723 + n\ + 6725 — n\ + Gng — ng — nsng) 

c = i{—22 + 12n5 — 2nl + 12n6 — 2ng — Sn^riQ + ngng + n^nl). 

Taking into account the admissible multiphcity of the resonances, we have 23 distinct 
cases of their positions. In 22 of those cases, however, the recursion relations for the 
coefficients of ( p.l| ) turn out to be inconsistent at the resonances (more details can 
be found in |^). The only good case is ni = 2,n^ = 2,125 = l^'^-e = 2, when 
n = -1,0,0,1,2,2,2,3,3,4,4,4, a = c = d = e = -3/2 and 6 = due to (U) 
and (U). In this case, we have to set f = g = —s/2 for the recursion relations to 
become consistent at n = 2, 3. Up to a scale transformation of variables, this is the case 
(|L7|) of the system (|1.6|). 



Three resonances at n = 0. In this case, resonances can lie in admissible positions 
only if go^o + ^0^0 = constant (we set constant = 1 w.l.g.) and 

a = /5 = 7 = 5 = — 1 a = —6 — b c = —6 — d — e {b — d){b — e) = 0. 

U b = d,e ^ d, we have 

{n + l)n^{n - 3)(n - A){n^ - Qn + 5 - 2d){n'^ - 6n + 5 - d - e)^ = 

i.e. six resonances lie in the positions n = —1,0,0,0,3,4. Denoting the positions of 
other six resonances as ni, 77,2, . . . , ?T-g, we obtain 

n2 = Q — 111 ris ^ rii 72,4 = 6 — 77,3 = = 

d = |(5 — 6ni + nl) e = 5 — d — 6213 + rig. 

There are four admissible choices of ni,n^. The cases rii = 2,713 = 1? ''^i = 3,ri3 = 1 
and rii = 3,n3 = 2 lead to inconsistent recursion relations for the coefficients of (|2.1| ). 
In the case rii = 1,723 = 2, when a = —6,6 = d = 0,c = e = —3, we have to 
set / = —2s,g = for the recursion relations to become consistent. Up to a scale 
transformation of variables, this is the case ( |1.8| ) of ( |l.(j| ). 
If 6 = e, we have 

{n + l)n^{n - 3){n - A){n^ - Qn + 5 - 2e){n'^ - Qn + 5 - d - ef = 
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i.e. six resonances lie in the positions n = —1,0,0,0,3,4. Denoting the positions of 
other six resonances as ni, ?t,2, . . . , ?t-6, we obtain 

n2 = Q — 111 = 6 — ris = n-^ = 



e 



|(5 — Grii + nf) d = 5 — e — Gn^ + n|. 



There are six admissible choices of rii, ^3. The cases ni = 2, 77,3 = 1, rii = 3, 71,3 = 1 and 
rii = 3,123 = 2 lead to inconsistent recursion relations for the coefficients of In the 

case rii = 1, 77,3 = 1, when a = c = —6, b = d = e = 0, we have to set f = g = —2s to 
make the recursion relations consistent at the resonances, thus obtaining the case ( |1.9| ) 
of ( |1.6|) . The case ni = 1, 713 = 2 with a = —6, b = e = 0,c = d = —3, where we have 
to set / = = —2s, leads us to the case ( |1.10|) of ( |1.6| ). The case rii = 2, 77,3 = 2 with 
a = —9/2,b = d = e = —3/2, c = —3, where we have to set / = 5^ = —s, leads us to the 
case ([TTTI) of (O). 



2.2.2. Case a + (3 = —2,7 + 6 > —2 We set go^o = 1 without loss of generality. If 
a 7^ /? or 7 7^ (5, there are two or more resonances in negative positions. Therefore we 
set a = P = —1, 7 = 5 > — 1 and obtain 

a = -6 - 6 d= -e + {2 + c)S ~3S^ + 

{n + l)n^{n - 3){n - 4:){n^ - 6n + 5 - 26)x (2.5) 

{n'^ -3{l-6)n + 2 + c-Q6 + 36^f = 

i.e. six resonances lie in the positions n = —1,0,0,0,3,4. Denoting the positions of 
other six resonances as ni,n2, . . . , uq, we find from ( p.5| ) that 77,2 = 6 — 77i, 771 = 1, 2, 3, 
b = |(5 — 672i +77^), 5 = 0, 773 = 775 = 1, n4 = 77e = 2, c = and d = —e. Then we check 
the consistency of recursion relations at the resonances. In the case ni = 1, we have to 
set e = g = and / = —2s at n = 1, and the system ( |1.6| ) becomes two independent 
Hirota equations 0. In the case rii = 2, we have to set e = (? = at n = 1 and / = — s 
at 77 = 3, thus obtaining two independent Sasa-Satsuma equations |]. In the case 
?7i = 3, the recursion relations are inconsistent at n = 3. 



2.3. Singular generic branches: h = 

When /i = 0, we set s = 1 without loss of generality. The consideration of dominant 
terms of ( |1.6| ) leads us to the following two cases to be studied separately: a + j3 = 
7 + 5 = — 1 and a + P = —1,7 + 6 > —1 (for the reason of symmetry, we omit 
a + P> -1,7 + 5 = -1). 

2.3.1. Case a + /3 = 7 + (5 = — 1 In this case, we have two possibihties: either two or 
three of the resonances lie in the position n = 0. 

Two resonances at n = 0. If both qqcjo and rgfo are some fixed nonzero constants, 
then the recursion relations turn out to be inconsistent at the resonance n = 1. 
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Three resonances at n = 0. We set go^o = i + ^fo^o w.l.g., and then find two 
possibilities: e = — 1 and e = 1. 
When e = —1, we have 



- 5 

-2-5a-5o2 h — ^ — -3«-3o2 



d = a — c 7 = a 



l+2a l+2« 

(n + l)n^{n - 2){n - 3)(n - 1 - k){n - 3 + k) = 

where k = —c — 2a, and only k = 0, 1, 2 are admissible. The case = 2 is related to 
k = through q ^ q,r ^ f,x ^ —x, t —t. In the case k = 1, the recursion relations 
are inconsistent at the resonances n = 2. In the case = 0, we have to set g = f at the 
resonance n = 1, and then the recursion relations are consistent at n = 2, 3 if and only 
if c = 2, c = —2, c = or c = 4, which, after scale transformations of variables, give us 
the cases (|1.12|) , (|1.13|) , (|1.14|) and ( |1.15|) of the system (|1.6|) , respectively. 
When e = 1, we have 

a = -^-5"-^°' b=-d= e = -a-c 7 = -! -a 

l+2a l+2a ' 

(n + l)n^{n - 2){n - 3)(n - 1 - l){n - 3 + /) = 

where I = c — 2a. In all the admissible cases Z = 0, 1, 2, however, the recursion relations 
are inconsistent at the resonance n = 2. 



2.3.2. Case a + (3 = —1, 7 + 5 > —1 We set goQ'o = i w.l.g., and then find the following 



two possibilities to have the resonances in admissible positions (see for more details): 

-2-5a-5o2 ^ _ _3a_3Q,2 



l+2a l+2a 
(l+2Q)e+(l+a)(7-7^)-a(i5-(5^) 



d 



(a7+(l+a)<5)e+7<5(2-7-5) 
(l+ci)7+a5 



(l+«)7+Q!5 

(n + - 2)(n - 3)(n + 7 + 5 + X 

(ra + 7 + 5- 2- ) = 



-2-5a-5o2 



(l+2Q!)e-(2ct+7)(5 ' 
(l+o)7+a5 

I, — 3q!— 3a^ 



1+2q: 1+2q 
.3/ 



d 



1 







(ra + l)n^{n - 2){n - ?>){n - 1 + c){n - 1 - c) = 0. 



:2.6) 



(2.7) 



In the case ( ^.6[ ), taking into account that four of the resonances should lie 
in positive integer positions and that 7, 5 should be integer, we have to set e = 
5 — 5^/ {l + 2a) and 7 = —5. Now the positions of resonances are n = —1, 0, 0, 0, 1, 1, 2, 3. 
At n = 1, we have to set g = 0,6 = —1 — 2a,a[a + 1) = 0, and then the recursion 
relations become consistent at n = 2, 3 as well. Both choices of a, a = and a = —\, 
lead to the case ( |1.16| ) of ( |1.6|) . 

In the case ( p.7|) , we have to set c = to ensure admissible positions of resonances. 
Then we obtain = at n = 1, and the system (|1.6|) becomes a pair of uncoupled 
equations. 
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2.4- Other branches 



We have proven that the system ( |1.6| ) admits good singular generic branches in the 
cases (|1.7| )- (|1.16|) only. In each of those cases, however, the system (|1.6| ) admits many 
other branches. They are Taylor expansions, which all are governed by the Cauchy- 
Kovalevskaya theorem in the case of ( [L.6D and need no analysis, as well as singular non- 
generic branches, which all have to be studied. All the singular non-generic branches of 
the cases ( [L.7|) - (|1.16|) of (|L6|) turn out to be good, in the sense that the exponents 
a, P, 7, S and the positions n of resonances are integer and the recursion relations 
are consistent. Omitting here the lengthy consideration of all cases, we give as an 
illustration the following two singular non-generic branches of the case (|1.7| ) (where 
h = l,a = -3/2): 

(i) a = —1, j3 = —1, 7 = —2, 6 = 2, qQ% = 4, positions of resonances are 
n = -4,-1,0,0,0,1,1,3,4,4,5,5; 

(ii) a = —2, /5 = 0, 7 = —3, 6 = 2, go^o = 8) positions of resonances are 
n = -5, -2, -1, 0, 0, 0, 2, 4, 5, 5, 6, 7. 



We can conclude now that the system ( |1.6| ) passes the Painleve test for integrability 
only in the cases ( |1.7| )-( P~!TBD . Let us establish in the succeeding section that the system 
]6|) in all the cases ([r7| ) -( [rTT6D possesses a Lax representation. 



3. Integrability 

3.1. Known cases and new cases 

First of all, let us note that, by the transformation 



X 

r 



t t 



f 



gexp(— icjj 
f exp(— ia;) 



in the cases 



■ X + pt 
gexp(ici;) 
r exp(iu;) 
ax + Tt 

with appropriately chosen real constants p, o, r, we can make s = f = g 
{^-{^J^) of (pD and / = ^ = in the cases ([LT2D-(|LT6D of (pj). 

If we set s = / = (? = in ( |1.8|) and ( |1.10|) , these two cases of (|1.6|) turn out to 
be simply related hj q —>■ q,q q,r f,f —>■ r. Since the case (|1.10|) is known to be 
integrable [^, the case ( |1.8|) is integrable as well. 

The integrability of the system ( |1.6| ) with ( |1.9|) , written in a form of coupled 
modified Korteweg-de Vries (mKdV) equations, was proven in (see also p7[|). 

A Lax pair for the case ( |1 . 1 1| ) of ( |1.6| ) was given in ^ . 

The integrability of the system (|1.6| ) with (|1.13|) was proven in [p9[1 . 

The cases ( |1.12| ), ( |1.14| ) and ( |1.15| ) are integrable due to the work 



m 



The remaining two cases, ( |1.7] ) and ( |1.16D , of the system ( |1.6| ) turn out to be new, 
and we prove their integrability in sections p.2| and |3.3|, respectively. 
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3.2. Case 

Without loss of generality, we set h = —1, a = —6, s = and study the case ( |1.7| ) of 
(|L6|) in the form of four coupled mKdV equations 

qt + Qxxx + Qqqqx + 6(grf)^. = 

qt + fcx + Qqqqx + 6(grf)3. = o 

n + J^xxx + 6rfr^ + 6(ggr)j, = 
n + ^xxx + 6rff^ + 6{qqf)^ = 0. 

The compatibility condition f/j — \4 + f/l^— Vf/ = of the linear problem = U"^, = 



V"^ becomes the system of two matrix mKdV equations |32 

Qt + Qxxx ~ <iQxRQ ~^ SQRQx = (3 2) 

Rt + Rxxx — SRxQR — 3RQRx = 



if we take the matrices U and V in the following block form 



Ho h [r 



V 




4/2 ; ^ 1^ 4i? ; V -2i?^. 2RQ 

I QxR — QRx —Qxx + 2QRQ \ 
\ —Rxx + 2RQR RxQ — RQx J 

where Ji and I2 are unit matrices, is a parameter. The choice 
\rq 




changes the system ( pl2|) into ( p.l|) . This proves that the new case ( |1.7|) of the coupled 
HONLS equations (|1.6| ) possesses a parametric Lax pair. The existence of a Lax pair 
indicates that the system is solvable under appropriate boundary conditions via 
the inverse scattering method. Indeed, we can solve the initial-value problem and obtain 
the iV-soliton solution of ( |3.1| ) by imposing appropriate constraints on the matrix- valued 
scattering data p6| . 

Usually, it is straightforward to obtain a multi-field extension of integrable two- 
component systems based on the Lax formulation. In fact, we can obtain an m- 
component generalization of the two-component system in the cases ( |1.8D , ( |1.1U| ) 



and ( |1.12| )-( [L15| ) by replacing two-component vectors in the Lax pair with m-component 



vectors. For the case ( |1 . 1 1[ ) , we have only to replace four-component vectors in the Lax 
pair with 2m-component vectors (this is clear if we express the Sasa-Satsuma equation 
in the form of coupled mKdV equations [§]). However, a multi-field generalization of 



the case (|1.7| ) as well as the case (|1.9|) is not straightforward. Thus we explain a multi- 



component generalization of the system (3J), which is analogous to that of the case 
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]9|) (cf |2^). For this purpose, we choose the matrices Q and R as 



Q 



Vol® I- Efc=r^ Vkl ®ek UqI®I- Efc=r^ UkCk ® I 



V vol0l~ El=i^ Vkl ®ek UOI0I + Efc=r^ UkCk ®I ) 

where / is the 2™^^ x 2™^^ unit matrix, and {ei, . . . , e2m-i} are 2™^^ x 2*""^ anti- 
commutative and anti-Hermitian matrices: 

Then the matrix mKdV equations ( p.2| ) become the system 

2m-l /2m-l \ 



fc=0 \ fc=0 

2 



/2m-i J = 0,1,..., 2m 



fc=0 V /c=0 / X 

Assuming that Uk and Vk are real and setting 

«2,-2 + i«2,-i = , = i,2,...,m 

W2j-2 + lf2i-l = 

we obtain a multi-component generalization of ( p.l| ): 

+ gj.xxx + Q^\(ik? qj,x + Q \ }^\rk? qA =0 

fc = l \fc = l /a; 

m / m \ 

fc = l \fc = l J X 

3.3. Case C^~TdO 

Without loss of generality, we set a = 2, s = 1, / = and study the case (|1.16| ) of (JlT 
in the form 



1,2, ... ,m. 



qt = \qxx + '2,[qqqx + rrqx + qrr^) ^^^^ 
n = ir^a: + 2(rfr^ + qqr^ + rqq^). 



Using the Mathematica package condens.m [^], we can check that the system ( |3.3| ) has 
two conservation laws for each rank from 1 to (at least) 4. The fact, that the conservation 
laws appear by pairs, is highly suggestive that equations ( |3.3| ) can be uncoupled by some 
transformation. Using the first conservation laws of ( ^.31 ), 

{qq)t = {i(«?xg - qqx) + {qqf + ^qqrf}^ 

{rf)t = {i{rxf — rfx) + (rf)^ + 2rfqq}x 
and introducing the new dependent variables by 

g = Mexp^i J vvdx'^ r = t>exp^i J uudx'^ (3.4) 
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we find that the system ( p.3| ) is equivalent to two independent Chen-Lee-Liu equations 

Ut = iUxx + 2uUUx Vt = iVxx + '2vvVx 

(here we have assumed that the dependent variables approach zero as x ^ Xq). Since 
a Lax pair for the Chen-Lee-Liu equation is well-known, we can straightforwardly 
obtain a Lax pair with two spectral parameters for the system ( |3.^j| ) by using a gauge 
transformation. This proves the integrability of the new case ( |1.16|) of (|1.6| ). It is 



remarkable that, due to ( [OD , the relations |g| = \u\, \r\ = \v\ are satisfied. These 
relations show that solitons observed in channel q and solitons observed in channel r 
interact only in their phase parts. 

4. Concluding remarks 

In this paper, we have carried out the singularity analysis of a system of two 
symmetrically coupled higher-order nonlinear Schrodinger (HONLS) equations with 
sufficient arbitrariness of the coefficients. Apparently, the form of coupled HONLS 
equations ( |1.6|) is more general than those which have been studied by means of the 
Painleve analysis |T^. Assuming the general form of coupled HONLS equations, 
we have exhaustively selected as many as ten cases in which the system passes the 
Painleve test for integrability well. The integrability of eight cases in the Lax sense 
has been already known, while the remaining two cases turned out to be new. We 
have verified the integrability of the two cases by considering a reduction of integrable 
matrix equations and a nonlinear transformation into uncoupled integrable equations, 
respectively. The work developed in this paper corresponds to the spirit of the work 
||T9| by Sofia Kovalevskaya, who originated the singularity analysis of nonlinear physical 
systems. 

We should comment on the result obtained in this paper. First, though the form 
of coupled HONLS equations (|L^) is general enough within limitations of computer's 
performance, we may obtain new integrable systems of coupled HONLS equations 
by assuming a more general form. Secondly, comparing the integrable cases (|1.2|)- 



( |L5|) of the single- component system with the integrable cases (L7)-( l.l(: ) of the two- 
component system, we notice that integrable one-component equations often have plural 
multi-component generalizations. We cannot identify each of ( |1.7|) -( prT6D with a multi- 



component generalization of only one of ( |L2|) - ([T75|) . For instance, the case (|1.7|) admits 
the reduction r — > to the Hirota case (|1.2|) as well as the reduction r —>■ q to the 



Sasa-Satsuma case (pTBI). Thirdly, study of exact solutions, which contain sufficiently 
many parameters, for integrable coupled HONLS equations has not been well-developed 
so far and deserve further investigation. 
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